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A great number of physical processes are described within the context of Hamiltonian scattering. 
Previous studies have rather been focused on trajectories starting outside invariant structures, since 
the ones starting inside are expected to stay trapped there forever. This is true though only for 
the deterministic case. We show however that, under finitely small random fluctuations of the 
held, trajectories starting inside Arnold-Kolmogorov-Moser (KAM) islands escape within finite time. 
The non-hyperbolic dynamics gains then hyperbolic characteristics due to the effect of the random 
perturbed field. As a consequence, trajectories which are started inside KAM curves escape with 
hyperbolic-like time decay distribution, and the fractal dimension of a set of particles that remain in 
the scattering region approaches that for hyperbolic systems. We show a universal quadratic power 
law relating the exponential decay to the amplitude of noise. We present a random walk model to 
relate this distribution to the amplitude of noise, and investigate this phenomena with a numerical 
study applying random maps. 

PACS numbers: 05.45.Ac 61.43.Hv 



Scattering is the general term used to describe sys- 
tems whose dynamics take place in an unbounded phase 
space, and such that the dynamics is trivial outside of a 
localised region called the interaction region. A classical 
example is a thrown particle which collides or interacts 
with a fixed target, then escapes from the neighbourhood 
of the obstacle. If the region where such particle inter- 
acts with the obstacle is much smaller than the whole 
space we locally assume that the particle comes from in- 
finity and is scattered again towards infinity. Alterna- 
tively, we may want to consider the dynamics of initial 
conditions starting inside the scattering region, in other 
words, the dynamics of an open system. We represent it 
by the dynamics, = f(x n ), taking place in an un- 

bounded phase space M such that, the scattering region 
is W C M. Accordingly, the orbits of initial conditions 
in W, the region of interest, may eventually leave W 
permanently. When the dynamics within a scattering 
region is chaotic, we say that we have a chaotic scat- 
tering. A typical signature for chaotic scattering is that 
characteristic quantities associated with the particles, or 
trajectories, after being scattered are sensitive to their 
initial conditions before entering the scattering region. 

Chaotic scattering is a very active research topic in 
dynamical systems. It has been found to describe a very 
broad class of dynamical processes. To mention a few 
cases, one may consider the dynamics of active chaos, 
i.e. chemical reagents or biological particles being ad- 
verted by the flow and at the same time undergoing in- 
trinsic transformations, such as biological reproduction 
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or chemical reactions. A remarkable example is found 
in the study of plankton population in the oceans ^ [2] ■ 
One may also consider the dynamics of blood flow and 
investigate the role of the chaotic advection in the resi- 
dence time of platelets and its consequence for the deposit 
and blocking of blood vessels [3] . A further example may 
be taken from the consideration of scattering particles 
in celestial dynamics [3], among many other dynamical 
phenomena. 

Because scattering dynamics are so closely related to 
scattering of particles in physical systems, we commonly 
refer to the dynamics of initial conditions in a region 
of the phase space as dynamical of particles whose dy- 
namics started in such region. From the dynamics point 
of view, chaotic scattering is characterised by the pres- 
ence of a chaotic non-attracting set containing periodic 
orbits of arbitrarily large periods as well as aperiodic or- 
bits distributed on a fractal geometric structure in the 
phase space — the chaotic saddle [5] . In analogy to more 
general dynamics, chaotic scattering can be hyperbolic 
and non-hyperbolic. In the former case, all orbits in the 
chaotic saddle are unstable, and a randomly chosen ini- 
tial condition has full probability of escaping within finite 
time. Hyperbolic scattering is also associated with an ex- 
ponential decay of probability for particles to be found in 
the scattering region after a given time t. That is, if we 
have a number N of randomly chosen initial conditions 
within the scattering region, the number of particles N(t) 
in that region decays exponentially on time and conse- 
quently 

P{t) (x e~ Kt , 

where P{t) = linijv(o)-»oo N(t)/N(0) is the probability 
for particles to remain in the interaction region after 
time t i and k a constant [5] . Non- hyperbolic scatter- 
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ing, on the other hand, is characterised by the presence 
of Kolmogorov-Arnold-Moser (KAM) islands in phase 
space. These islands surround marginally stable periodic 
orbits. Orbits from the outside can spend a long time in 
the vicinity of the KAM islands before escaping, and this 
"stickiness" effect causes a slower escape dynamics than 
that found in hyperbolic systems: non-hyperbolic escape 
is characterised by a power law distribution of probabil- 
ity, 

P(t) oc r a , 

for large enough t [6-8 . An important characteristic is 
that due to the invariance of the islands and the area 
preserving property of the dynamics, if initial conditions 
are chosen inside an island, orbits are expected to be 
trapped there forever [5] . 

The above discussion refers to completely determinis- 
tic dynamics. However, since most phenomena in nature 
are always subjected to small fluctuations from their sur- 
rounding environment, a central question is: how does 
the dynamics of scattering systems change in the 
presence of random perturbations? 

Here we focus on random perturbations of the param- 
eters defining the system. From the Dynamical Systems 
perspective, this situation can be formalised by the con- 
cept of random maps, which we will apply throughout 
this paper, 

where we randomly choose slightly different maps /„ for 
each iteration n (see Eq. [2] below) . It is known that such 
dynamics has well-defined (in the ensemble sense) values 
of dynamical invariants such as fractal dimensions and 
Lyapunov exponents O HH [17] . Note that we associate 
the choice of the map with the iteration. Therefore, all 
initial conditions in a given iteration are mapped by the 
same sequence of random maps. On the other hand, if 
we implement independent random noise added to each 
trajectory, what has been previously considered for the 
random noisy dynamics |10H14j , the situation is very dif- 
ferent, as we shall discuss later. 

Let us recall that hyperbolicity is a structurally stable 
characteristic of dynamical systems, i.e. hyperbolic sys- 
tems are robust under small smooth perturbations of the 
system. But on general grounds, one expects the pertur- 
bation of a non-hyperbolic dynamical system to result 
in a hyperbolic dynamics. We remark that the notion 
of structural stability is related to perturbations of the 
whole dynamical system and not to noisy diffusion-like 
perturbations of individual trajectories [T5] . Although it 
is possible to define random invariants in more general 
grounds [5], for our context it makes no sense to talk 
about dynamical invariants such as the Lyapunov expo- 
nent and fractal dimensions in the noisy dynamics in this 
latter case, because trajectories will be "smoothed-out" 
in small scales and all fine-structure dynamical structures 
will disappear. In the case of perturbation of parameters, 




FIG. 1: (Colour online) Simplified representation of the orbits 
on the torus. In (a), it is represented a single torus and its 
centre. In (b), it is represented the effect of the perturbation. 
The difference between the dashed and the continuous line is 
due to the effect of the noise, that deslocate the centre of the 
torus from 0\ to O2. In (c) we show what would be expected 
for 4 iterations. Each continuous line represents the orbit of a 
particle on the torus for each iteration, hence different values 
of perturbations at each iteration. The centre is expected to 
move around Oi at each iteration and for long enough time, 
the distribution of centres would fill densely the area within 
the dashed circle. 

on the other hand, we have "random maps" [16], which 
are known to have well-defined dynamical invariants in a 
measure-theoretical sense [17], 

A further observation is that, since hyperbolic systems 
are structurally stable, we expect that small perturba- 
tions would not result in qualitative changes to the dy- 
namics of hyperbolic maps. In particular, escape should 
continue to be exponential. For non-hyperbolic systems, 
the same can not be expected. Thus, we expect to see 
qualitative changes due to the perturbations. As an ex- 
ample, we could heuristically see the KAM tori no longer 
as confining, since now the perturbations can always push 
orbits out of the region of an island, and conversely or- 
bits on the outside can be pushed inside a torus by the 
perturbations. The destruction of this essential feature 
of non-hyperbolic scattering leads us to hypothesise that 
even arbitrarily small random perturbations to the dy- 
namics will turn a non-hyperbolic scattering system into 
a hyperbolic-like dynamics. Hence, we expect the prob- 
ability distribution of escaping time for perturbed non- 
hyperbolic systems to become exponential. 

In this paper we explore and test this hypothesis us- 
ing a simple scattering map as a model system. We de- 
rive an analytical theory to predict the scaling of the 
escape rate of non-hyperbolic scattering systems under 
weak perturbations. The theory predicts that the (expo- 
nential) escape rate goes to zero as the amplitude of the 
perturbation decreases following a quadratic law, which 
we argue is universal for all perturbed non-hyperbolic 
systems. We verify this prediction numerically, and find 
that it describes well the results of our simulations. 

In order to better understand the effects of random pa- 
rameter perturbations on the dynamics of non-hyperbolic 
Hamiltonian systems, we develop now a simple statistical 
model of the motion of a single trajectory in the presence 
of perturbations. 

In the absence of random perturbations, the islands act 
as barriers in the phase-space; Due to the area-preserving 
property of such dynamics, those invariant structures can 
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FIG. 2: Phase space for the map [2] for A = 6.0. The inset 
shows a blowup of the region (x, y) € [2.05, 2.20] x [0.44, 0.47]. 



not be transposed [5] , thus trajectories trapped inside is- 
lands never escape. Now suppose a small noise is added 
to the system's parameter. The effect of noise is to vary 
randomly the parameter from iteration to iteration, but 
noise amplitude is assumed to be small enough so that 
the parameter lies in the range where the system is non- 
hyperbolic. The effect of one iteration could be, say, to 
change the parameter A to a slightly smaller value, but 
still keep it within the range such that the dynamics is 
non-hyperbolic. The effect of the new parameter is then 
to cause a small shift in the KAM structures compared 
to the value of A in the previous iteration. For the next 
iteration another value of A is chosen, and the KAM is- 
lands again move slightly in phase-space. This happens 
at every iteration, and the global effect is a sort of ran- 
dom walk of the KAM structures around their positions 
in the absence of noise. Fig [T] is an idealised illustration 
of the orbit of a particle on a torus. We represented three 
steps, a simple torus - Fig [I] (a)-, the effect of the new 
parameter - Fig[l](b) -, and finally, what would be ideally 
expected for 4 iterations - Fig [T] (c) . 

In other words, the perturbations can be imagined to 
cause orbits to gain motion in the direction transversal to 
the tori. The magnitude of this transversal component of 
the motion is proportional to the intensity of the pertur- 
bation (in the case of small perturbations) . Since only the 
transversal component can cause an orbit started within 
a KAM island to escape, we focus on this component of 
the motion alone. 

The transversal motion is due entirely to the random 
perturbations, and as a first approximation we model it 
as a one-dimensional random walk. The size of the step 
5 of the random walk is proportional to the amplitude 
of the perturbation. After n steps, the typical distance 
D from the starting position reached by the walker is 
D ~ y/n8 [18] . Let Do be a typical transversal distance a 
particle needs to traverse in order to cross the last KAM 
surface and escape. Thus the average time (number of 



steps) r it takes for particles to escape scale is given by 
Dq ~ tS 2 . So t scales as t ~ 5~ 2 . The conclusion is 
that our simple model predicts an exponential decay of 
particles, with a decay rate k ~ r _1 scaling as 

k~6 2 . (1) 

Note that the dynamics described by random maps 
here is very different from that where random noise is 
added to individual trajectories. In the latter case, an 
analysis in terms of stochastic stability is more adequate, 
because the behaviour of, for example, two distinct or- 
bits evolved from different initial conditions cannot be 
compared in terms of (Lyapunov) stability. In fact, even 
by following two distinct time evolutions from the same 
initial condition we may observe completely distinct be- 
haviour due to the two different noise realisations, pos- 
sibly not correlated. In particular, for scattering under 
the presence of noise [T2] , it has been shown that the ex- 
ponential escape rate is preserved in weakly dissipative 
systems and that also an exponential escape rate should 
be expected for non dissipative ones, for noise intensity 
above a threshold. Indeed, because trajectories under 
different noise realisations may not be correlated for a 
certain level of noise, we expect a resultant random be- 
haviour to be dominant, and therefore a distribution of 
random values of escaping time. Hence, an exponential 
distribution naturally follows. This scenario has already 
been pictured in [12] by considering that above the crit- 
ical level of noise, an effective random redistribution of 
particles results in a blurring of the fine structure of the 
KAM tori. In the dynamics given by random maps con- 
sidered here, on the other hand, all the trajectories are 
evolved under the same realisation and fine-scale struc- 
tures are preserved for each map, allowing a well-defined 
analysis of structural properties of limit sets [3] ■ In this 
case, one can, for example, evaluate the Lyapunov stabil- 
ity of orbits by comparing the convergence towards one 
another or the divergence of different random orbits gen- 
erated by the evolution of distinct initial conditions, since 
they are subject to the same sequence of perturbations. 
It is also possible to define consistently the fractal dimen- 
sion of invariant sets under the iteration of the random 
maps [17], unlike in the case of independent realizations 
of noise for different particles. Intuitively, we can think of 
such a system as perturbations of the dynamics itself, for 




FIG. 3: KAM structures under the perturbed dynamics for 
different values of f. Fig^a), £ = 0, Fig^b), £ = 0.0004, 
and Figgc), £ = 0.002. 
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FIG. 4: (colour online) Probability distribution of escaping time t from the region W = {\x\ < 5.0, and \y\ < 5.0} for 
different values of £. Initial conditions were randomly chosen in the line x £ [2.05, 2.07], y = 0.465 inside nest KAM structures. 
For each value of f, we present the exponent 7 that best fits the exponential region of the probability distribution. Fig[4|a) 
shows all used values of £, FigHfb), shows small values of t. 



example, a set of particles that evolve subject to a field 
that undergoes random fluctuations. Therefore, all par- 
ticles are thought to evolve under the same sequence of 
dynamical laws, rather than each particle being subject 
to random uncorrelated kicks. This key characteristic is 
what allows us to consider our model based on the ran- 
dom walk of the tori and develop our approach in order 
to derive the power law described by equation [T] As we 
shall see in what follows, it also allows us to consider 
singular sets of initial conditions, which we will use to 
compute the fractal dimension. Furthermore, we remark 
again that the initial conditions used here are always cho- 
sen inside the original KAM structures. For this set of 
particles, if no random perturbation is considered, there 
should not have any escape at all due to the are preserv- 
ing property of our map. Simply the fact that such escape 
indeed happens is already by itself a novelty characteris- 
tic introduced by the perturbations of the dynamics and 
explained by our theory. 

In order to test our theory, we numerically obtained the 
probability distribution P(t) that particles have not es- 
caped at time t from a given region XV of the phase-space 
under the dynamics of an area-preserving non-hyperbolic 
map. We have chosen the map 

x n +i = K[%n ~ (x n + y n ) 2 /4] 

1 (2) 
Vn+l = + {%n + ynf/4], 

An 

which is non- hyperbolic for A < 6.5 [7]. We define 
W = {\x\ < 5.0, \y\ < 5.0}. The initial conditions 
were randomly chosen with uniform probability in the 
line x £ [2.05, 2.07], y = 0.465. For this interval, the 
particles start their trajectories inside a KAM structure 
(See Fig [2]). Because we are interested in perturbing the 
system, we chose A* = 6.0. Then, for each iteration n 
we randomly chose a perturbation |e n | in the interval 
|e n | < £, where £ is the amplitude of the perturbation. 



Indeed, for the simulations we observe, as expected, the 
scenario pictured in Fig. [T] We show in Fig. [3] a portion 
of the phase space containing KAM islands for different 
values of perturbations. 

We show in Fig [4] the probability distribution P(t) for 
different values of £. We clearly see that the escape fol- 
lows an exponential law in a broad range of times, for 
various perturbation amplitudes. For each value of £, we 
show the exponent 7 that best fits the exponential law 
Pit) ~ e~ 7 * in the range of times for which the proba- 
bility distribution is exponential. 

We also observe that for small values of the perturba- 
tion amplitude, the stickiness plays an important role on 
the escaping time distribution of particles that leave the 
non- attracting chaotic set. This is because even after a 
particle escapes from a KAM island, it may be trapped 
for long times by Cantori. Under small amplitude of 
perturbations, this causes a slower escape than that pre- 
dicted by our model. This signature of the stickiness is 
observed in the tail of the distribution shown in Fig. [4] as 
a cutoff for long times in the exponential escape of par- 
ticles. The exponential probability distribution of escape 
time agrees with our hypothesis that the dynamics of the 
perturbed system is hyperbolic-like in the presence of the 
random perturbations. 

Using our random walk escaping model, from Eq. [T] 
and Fig. [5j we identify k = 7, and 8 = £. In order to 
construct our model, we have not taken into account the 
effect of stickiness. Therefore, it is reasonable to expect 
that the obtained exponent should actually be smaller 
than 2. We show in Fig. [5] the dependence of 7 on the 
amplitude of the perturbation £, for the map given by 
Eq. 1(2). First of all, we notice that 7 indeed follows a 
power-law, as our model predicts. We also obtained a 
good agreement with the square law predicted by Eq. [l] 

Another fundamental characteristic, differentiating 
non-hyperbolic chaotic scattering from the hyperbolic 
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case, is the fractal dimension of singular sets. One ex- 
pects that choosing initial conditions from a line in W, 
the set of particles that remain in W after a given time 
To should form a Cantor set with fractal dimension d < 1 
in the case of hyperbolic dynamics. It is due to the expo- 
nential decay of probability characteristic of hyperbolic 
systems. On the other hand, the algebraic decay of non- 
hyperbolic scattering leads to the maximal value of the 
fractal dimension, d — 1 7J. However, it is well known 
that in practice one would need in most cases to go down 
to unreasonably small scales to sec the actual d = 1 pre- 
diction in a non-hyperbolic scattering system. For scales 
that are physically meaningful, in most cases one finds 
that the chaotic saddle and its associated sets (stable 
and unstable manifolds, etc.) have a position-dependent 
effective dimension 2, 8J, which is lower than 1. 

As an example of scattering function, we consider the 
time-delay T(x). The idea is to calculate the fraction 
/(e) of pairs of initial conditions started in random po- 
sitions of the phase space with distance e of each other 
whose trajectories diverge within the scattering region, 
and whose orbits end up escaping through completely dif- 
ferent routes. In order to compute that, we chose initial 
conditions inside a KAM structure. We chose yo — 0.465, 
and different values of xo were randomly chosen to be- 
long to the interval [2.05,2.15]. For a fixed value of un- 
certainty e we chose xq and compute \T(xq) — T(xq + e)\. 
If \T(xq) — T(xq + e)| > 2 we say that x$ is e-uncertain. 
Dividing the number of e-uncertain points by the total 
number of initial conditions we obtain /(e), the uncer- 
tain fraction. It can be shown that for chaotic scattering 
systems /(e) scales as a power law e' 3 , and that the box- 
counting fractal dimension d is given by d = 1 — /3 |19j . 
We expect that the smaller the amplitude of noise, the 
closer the d is to 1, which would correspond to the non- 
hyperbolic limit. We show in Fig.[6]the estimated fractal 
dimension of the function T(x), for different values of £. 




FIG. 5: Different values of exponent 7 as a function of £ from 
Fig. [4] as well as the power law that best fits the distribution. 
In the inset, we consider only values of £ < 0.02. 



□ ^ = 0.1024 
o E, = 0.0512 




FIG. 6: (colour online) Estimated fractal dimension of the T, 
for different values of £. We notice that, when the amplitude 
of the perturbation is decreased, the dynamics approach the 
non- hyperbolic limit, and so does the estimated values of frac- 
tal dimension. The inset shows the value of (3 as a function 

ofe 



We notice that when the amplitude of the perturbation 
is decreased, the dynamics approach the non-hyperbolic 
limit, and so do the estimated values of fractal dimension, 
which approach 1 as £ — > 0. Furthermore, we notice that 
/? scales as a power law with the amplitude of noise, i.e. 

p ~ c ): - 7 . 

In conclusion, we have shown that random small per- 
turbation of non-hyperbolic maps leads to hyperbolic be- 
haviour. We have presented a random walk model to ex- 
plain the escape of particles from inside invariant KAM 
islands. Such particles are expected to be trapped there 
forever in the absence of perturbation. When the per- 
turbation is present, we show that, not only particles are 
able to escape, but also they have hyperbolic behaviour, 
with a universal quadratic power law relating the expo- 
nential decay to the amplitude of the perturbation. We 
further investigated the hyperbolic behaviour estimating 
the fractal dimension of the set of particles that remain 
inside the islands for a given time. We show that indeed 
the fractal dimension is in agreement to what is expected 
for a hyperbolic system. 

We call the attention of the reader to an important 
issue on scattering dynamics clarified by our results. 
Note that despite the generality of scattering phenomena, 
mostly of the modelling of real processes has assumed a 
simpler chaotic dynamics, that is, given by hyperbolic 
structures. Nevertheless, it is widely argued that most of 
the flows experimentally significant are non- hyperbolic; 
see references therein [2]. In active dynamics, or blood 
flow, for example, the modelling describe systems where 
we have trapped vortices. Therefore, in the absence of 
perturbations the particles would not leave the trapping 
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regions; This would correspond to plankton not being 
able to leave a given region in the sea, or blood particles 
being trapped indefinitely in vortices, which is in contrast 
with most observations. As a consequence, although the 
dynamics is described by non-hyperbolic scattering, the 
escape is in many cases assumed to follow an exponen- 
tial law, without any further explanation [TJ [3J. Since 
it is reasonable to assume that these flows naturally ex- 
perience random perturbations, our theory provides an 
important bridge to allow one to encompass natural phe- 
nomena under the non-hyperbolic scattering framework. 

Apart from striking implications to the above men- 
tioned dynamics, our result may play a fundamental role 
on the understanding of how sources, and escaping rates 



of particles act on the dynamics of Saturn Rings [3] . The 
rings are thought to be marginally stable periodic orbits, 
and the gaps to be rational tori or made up by commen- 
surable frequencies. They are subjected to gravitational 
perturbation and perturbations of the electromagnetic 
field, due to solar storms, solar winds, magnetic storms, 
etc. [3] , or due to the dynamic of ring current around Sat- 
urn [3j . Such random perturbations would be represented 
by our perturbation on the parameter A in Eq. [2] Hence, 
if measured, we would expect an exponential-like escape 
of particles. Therefore our results may shed some light 
onto the understanding of the escape of small particles 
like in the Saturn's E ring [3]. 
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